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Abstract 

We study the effect of the induced interaction on the superfluid transition temperature of a 
Fermi gas with a BEC-BCS crossover. The Gorkov-Mehk-Barkhudarov theory about the induced 
interaction is extended from the BCS side to the entire crossover, and the pairing fluctuation 
is treated in the approach by Nozieres and Schmitt-Rink. At unitarity, the induced interaction 
reduces the transition temperature by about twenty percent. In the BCS hmit, the transition 
temperature is reduced by a factor about 2.22, as found by Gorkov and Mehk-Barkhudarov. Our 
resuh shows that the effect of the induced interaction is important both on the BCS side and in 
the unitary region. 
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I. INTRODUCTION 



One of the most important developments in experiments on ultra-cold atoms is the ob 



servation of BEC-BCS crossover 
strongly-coupled superconductors 



which was originally predicated for 



lOj. In experiments on ultra-cold atoms, the interac- 
tion between atoms can be tuned by the technique of Feshbach resonance, and the system can 
evolve smoothly from a Bardeen-Cooper-Schrieffer (BCS) pairing state to a Bose-Einstein 
condensation (BEG) state of diatomic molecules. In the BCS limit where the interaction 
is weakly attractive, the atoms are paired into a BCS state below a critical temperature, 
very similar to electrons in conventional superconductors. In the BEC limit where the in- 
teraction is weakly repulsive, tightly-bound diatomic molecules are formed, and a molecular 
.3EC state appears below a critical temperature. Near the resonance, in the unitary region 
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12| where the size of the scattering length is much larger than the inter-particle spacing. 



the system is strongly correlated in both the normal and superfiuid states. 



The BEC-BCS crossover can be qualitatively explained by a mean-field BCS theory 10|. 
In this theory, the size of atom pairs decreases as the system goes from the BCS side to 
the BEC side. In the BEC limit, the pair size is so small that atom pairs become diatomic 



molecules. However, the mean-field theory 



transition temperature in the BEC limit 



3redicts an exponentially-divergent superfiuid 



131 , which is against the result from the theory 



about an ideal Bose gas. Nozieres and Schmitt-Rink (NSR) IJ] first pointed out that the 
pairing fluctuation must be taken into account to obtain the correct superfiuid transition 
temperature of the BEC-BCS crossover. The pairing fluctuation is especially important 
in the BEC limit where nearly all atoms become thermal molecules at T^. Following NSR's 
pioneer work, many theoretical studies have focused on irnproving NSR's method and ex- 



tending their analysis to the broken symmetry state 
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which was recently reviewed in Ref. 
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2^, 



26|. 



However in the BCS limit, a different type of fluctuation is important. Gorkov and Melik- 
Barkhudarov (GMB) found that there is a modification to the pairing interaction due to 



the many-body medium, referred to as the induced interaction [28| . The induced interaction 
suppresses pairing considerably and reduces the superfiuid transition temperature Tc by 
a factor about 2.22 with respect to the mean-field in the BCS limit. The fluctuation 
considered by GMB is in the particle-hole channel, different from that in the particle-particle 
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FIG. 1: Critical temperatures of BEC-BCS crossover. The solid line is our result after taking into 
account the induced interaction in the NSR approach, the dashed line is the original NSR result, 
the dotted line is the result from BCS mean-field theory, and the dot-dashed line is GMB's result 
given by Eq. These theoretical results are compared with the data from QMC simulations 



31 



32l ] shown in symbols. 



channel considered by NSR. In the BCS limit the NSR fluctuation is much less important 
than the GMB fluctuation. When the system moves from the BCS side towards the BEC 
side, the GMB fluctuation becomes weaker and the NSR fluctuation becomes stronger. In 
the BEC limit, the NSR fluctuation is dominant. In an accurate description of the BEC-BCS 
crossover, both GMB and NSR fluctuations should be treated pro per ly, which has not been 
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addressed except in two recent renormalization-group studies 

In this work, we present our result about the induced interaction in the whole BEC-BCS 
crossover. First we extend the GMB theory about the induced interaction from the BCS 
limit to the strongly-interacting region. Then we consider the induced interaction in the NSR 
framework and compute T^. for the entire BEC-BCS crossover. Our main result is shown 
in Fig. [H Compared with the original NSR result, the superfluid transition temperature 
Tc is reduced considerably both on the BCS side and in the unitary region. In the BCS 
limit, we recover the GMB result. At unitarity, the critical temperature T^. is found to be 



Tc = 0.178Tp, close to the results from Quantum Monte Carlo (QMC) simulations (3ll. l32|. 
and about 20% smaller than the NSR result. Our result shows that induced interaction plays 
an important role in the unitary region and on the BCS side. Discussions and conclusions 
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are given in the end. 



II. THE INDUCED INTERACTION 

A Fermi gas with a wide Feshbach resonance can be described by a single-channel model, 

where the coupling constant is given hj g = Anh'^as/m, is the scattering length, and ipa- 
is the field operator for spin component cr. In this work, we consider only the homogeneous 
spin-balanced case where the densities of spin-| and spin-| atoms are the same. 

In the BCS limit when the interaction is weakly attractive, Gorkov and Melik- 
Barkhudarov (GMB) 271] showed that in the particle-hole channel there is a correction 
to the pairing interaction from the many-body background, given by the Feynman diagram 
shown in Fig. [2](a). Other diagrams of the same order are not as important. For example, 
the diagram in Fig. [2](b) corresponds to an effective interaction between atoms with the 
same spin component, which is strongly suppressed at low temperatures. 

Beyond the BCS limit, when the interaction is strong, higher order diagrams are impor- 
tant. We generalize the GMB approximation by considering all the diagrams of the same 
general type as the GMB diagram shown in Fig. [2](a) and summing all these diagrams to- 
gether as shown in Fig. W^c), and obtain the induced interaction in the normal state given 
by 

Uind{Pl,P2;P3,P4) = -— 7 T, (2) 

1 + 9X{Pi - Pi) 

where pi = (kj,u;;J is a vector in the space of wave-vector and frequency, and ui = {21 + 
1)11 /{h(3) is the Matsubara frequency of a fermion, /3 = 1/ {ksT). The total momentum and 
energy are conserved in the scattering, Pi + P2 = + P4- The function x is taken as the 
polarization function of a non-interacting Fermi gas with the same chemical potential /i, 
given by 



X{p') = ^J2^oip)Qoip + p') 

p 

d^k /k - /k+k' 



(27r)3 ih^i + Ck - ek+k' 



(3) 
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FIG. 2: Feynman diagrams in the particle-hole channel, (a) The induced interaction considered by 
GMB. (b) An effective interaction between atoms with the same spin component, (c) The induced 
interaction in the generalized GMB approximation. 

where = (k', Qi = 2111 / {h(3) is the Matsubara frequency of a boson, V is the volume, 
/k = + exp(/3^k)] is the Fermi distribution function, = Ck — /W, and ek = h'^k'^/2m. 
The Green's function of a non-interaction Fermi gas Qo{p) is given by Qo{p) = h/{ifvjJi — ■Ck)- 
Including the induced interaction, the effective interaction between two atoms with dif- 
ferent spin components is given by 

UtotiPuP2;P3,P4) = 9 + Umdipi,P2;P3,P4) 

= ^ . (4) 

1 + 9X{Pi - Pi) 

Although the effective interaction is a function of transferred momentum and frequency, at 
low temperatures only its s-wave part plays an important role on pairing. As in GMB's 
work, we approximate this s-wave component g' by averaging the polarization function 

When yU > 0, the average of the polarization function (x) is obtained by setting the frequen- 
cies and total momentum to zero and taking all the initial and final states of atoms from 
the Fermi surface, i.e. ki = — k2, ks = — k4, ki = k2 = = k^^ = kp, which yields 



m 



f-l POO ^ 


k' 


-2k 


1 dcos6 J dk—fk\n 


k' 


+ 2k 



(6) 



where k' = |ki — k4| = kp^J2{l + cos6'), kp is the Fermi wavevector, and 6 is the angle 
between ki and When the chemical potential turns negative on the BEG side, the 
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Fermi surface disappears, and the average (x) is taken at zero frequency and in the hmit 
that all the momentum go to zero, ki = k2 = = 0, same as /c' limit of Eq. ([6]), 

{x) = -^rdkf,. (7) 



2n^h^ JO 

In both cases, the function (x) is always negative and monotonically decreasing with the 
increase in the chemical potential fi. 

In the BCS limit, the critical temperature is much less than the Fermi temperature 
Tp. Near Tc, one obtains 

where A/'(ei?) = mA;i?/(27r^^^) is the density of states for one spin species at Fermi energy. 
The effective s-wave interaction g' is approximately given by 

With the effective pairing interaction g', the GMB result of can be obtained, 

/ 9 \ 7/3 

yJGMB) ^ / f \ J_j.^^n/2k^a. ^ 0.28Tire"/2fc^"% (9) 

where 7 = e*^, c is the Euler constant. The GMB result j'^^*^^^) jg smaller by a factor of 
(4e)^/3 ^ 2.22 than the mean-field T^. 

III. THE T-MATRIX AND CORRECTION TO DENSITY 

To obtain the critical temperature for the whole BEC-BCS crossover, we compute the 
T-matrix, as shown in Fig. [3], 

t(p') = = (10) 

l+9'Xp{p') l/9+{x) + Xp{p') ^ ' 

where the pair susceptibility Xp{p') iii particle-particle channel is given by 

Xpip') = -^^^Go{p)Go{p' - p) 

p 

d'k h + /k'-k - 1 . . 

(27r)3zni],-ek-^k'-k' ^ ^ 

Comparing with the conventional T-matrix approach, we have replaced the coupling constant 

g by the effective s-wave interaction g' due to the induced interaction. 
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FIG. 3: Diagrams of the T-matrix. The wiggled line represents the effective s-wave interaction g' . 

According to Thouless criterion, the superfluid instability at is due to the divergence 
of t{p' = 0) which is equivalent to 

m , , f d^k /2/k-l 1 \ 

■ ix)= 77^A + 7T- h (12) 



4^2™, J (27r)3 V 2^k 2ek^ 

where the last term on the right-hand side is the counter term due to vacuum renormaliza- 
tion. Comparing the T^-equation given by Eq. f|T2l) with that in the BCS mean-field theory, 
the effect of the induced interaction is equivalent to making the scattering length larger. In 
the BCS limit, the induced interaction leads to a reduction of Tc from the mean-field result, 
as given in Eq. In the EEC limit, the effect of the induced interaction is negligible. 

As Nozieres and Schmitt-Rink pointed out [l^ , pairing fluctuations in the particle-particle 
channel are important especially on the BEC side. In the NSR theory, the total atom density 
includes not only the fermion density in the mean-field approximation, but also contributions 
from fluctuations of molecular fields. In the BEC limit, near Tc, the fluctuation contribution 
is dominant and the critical temperature Tc is given by the BEC temperature of an ideal 
Bose gas. One way to take into account the NSR effect is to calculate the Hartree self-energy 



generated by the T-matrix and its contribution to density [2l| . With our T-matrix given by 
Eq. (ITO!) which includes the induced interaction, the self energy is given by 

^(p) = ^Y.^(p'w-p). (13) 

p' 

To the first order, the Dyson's equation is given by 



and the particle density is given by 

2 



n 



h(5V 

where the mean-field density is given by 



J]6?(p)e-'"'°^ =n/ + An, (14) 



■(27r)= 



and the fluctuation contribution An is given by 

2 



An 



(15) 



If we omit the (x) term due to the induced interaction in the T-matrix given by Eq. (fTOll . 
the density equation given by Eq. f|T5|) is the same as that in the NSR theory. 



IV. THE SUPERFLUID TRANSITION TEMPERATURE 

The superfluid transition temperature Tc as a function of the total atom density n can 
be solved from the two coupled equations f|T2|) and f|T^ . In the BCS limit, since rif ^ An, 
the GMB result about Tc given by Eq. can be recovered. In the BEG limit, at Tc, the 
mean-fleld density is negligible. An ^ n/, the T-matrix is proportional to the propagator of 
noninteracting molecules, and the density of total atoms is approximately given by density 
of molecules. Thus the transition temperature in the BEG limit is given by the condensation 
temperature of an ideal Bose gas, = 0.218Tp. 

Our numeric result of the critical temperature for entire crossover is shown in Fig. [H 
As expected, in the BGS limit, it agrees with GMB theory; in the BEG limit, it recovers 
the condensation temperature of ideal molecules. At unitarity, we obtain Tc = O.lTSTp, 
which is close to the QMG result Tc = 0.15(l)Tp 3l|, |3^. In comparison, the results from 
other theoretical studies are Tc = 0.222Tp in the original NSR theory, Tc = 0.160Tpin 
a full self-consistent NSR treatment |23|, Tc = 0.26Tp in pseudogap crossover theory 2l|, 



Tc = 0.264Tp 29| and Tc = 0.13Tp 30|] in renormalization group studies. Gompared with the 
original NSR result, our critical temperature is about 20% lower, implying that the induced 
interaction still plays an important role in the unitary region. The chemical potential at 
Tc in our results is /x(Tc) = 0.598Tj., higher than QMG results, fi{T^) = 0.493(14) [Sll and 



0.43('ll 



|. Our results can probably be improved by self-consistently taking into 
account the self energy in the computations of the induced interaction and T-matrix. This 
issue will be addressed in our further studies. 
As reported in previous works 14 , 1^, 



20 



231], we also flnd that the critical temperature 



reaches a maximum on the BEG side, as shown in Fig. [TJ Gompared with the original NSR 
result, the position of this peak is further away from the resonance due to the induced 
interaction. Our results show that the peak is located at l/kpCLs = 0.437, and Tp^*^^ = 
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0.231Tf, close to the QMC estimation of the peak position l/kpas > 0.474(8) and TP""^ > 
0.252(15)Tf [311. 

The effect of the fluctuation in the particle-hole channel on the superfluid transition 
temperature was also studied in the renormalization group approach for a two-channel model 
mostly in the wide resonance case 29|] and for a single- channel model 30|. In the BCS limit, 
the superfluid transition temperature was found in agreement with the GMB result in Ref. 



291] , and smaller than the GMB result in Ref. 30| due to the simplification in the momentum 



dependence of the interaction vertex. At unitarity, the superfluid transition temperature 
was found to be = 0.264Tf ^ and 7; = O.lSTp [aoj, while we obtain = 0.178Tp 



and the QMC result is = 0.15(l)rp |3l|, On the BEG side when kpa^ = 0.5, the 



291], and our result shows 



superfluid transition temperature was found to be ~ 0.25Tp 
Tc ~ 0.22Tp. These quantitative differences may be resolved in future studies with better 



theoretical treatments. 



V. CONCLUSION 



In conclusion, the effect of the induced interaction due to the many-body medium is 
studied in a Fermi gas with the BEG-BGS crossover. The GMB theory is extended from the 
BGS limit to the entire crossover. With the induced interaction considered, the superfluid 
transition temperature Tc is computed for the entire crossover in the NSR framework. The 
induced interaction reduces the critical temperature considerably on the BGS side and in 
the unitary region. Our results of = O.lTSTp at unitarity and the T^-peak location are 
in reasonable agreements with results from quantum Monte Garlo simulations. Our results 
show that the effect of the induced interaction is important both in the unitary region and 
on the BGS side. We would like to thank T.-L. Ho for helpful discussions. This work is 
supported by NSFG under Grant No. 10674007, and by Ghinese MOST under grant number 
2006GB921402. 
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